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Abstract 

The equations for the second-order gravitational perturbations produced by a compact-object 
have highly singular source terms at the point particle limit. At this limit the standard retarded 
solutions to these equations are ill-defined. Here we construct well-defined and physically mean- 
ingful solutions to these equations. These solutions are important for practical calculations: the 
planned gravitational-wave detector LISA requires preparation of waveform templates for the ex- 
pected gravitational-waves. Construction of templates with desired accuracy for extreme mass 
ratio binaries, in which a compact-object inspirals towards a supermassive black-hole, requires 
calculation of the second-order gravitational perturbations produced by the compact-object. 
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In an extreme mass-ratio binary, where a compact-object (CO) of mass \i (e.g. a neutron 
star or a small black hole) orbits a supermassive black hole of mass M, the interaction 
between the CO and the gravitational perturbations that it produces exerts a force on 
the CO. This phenomenon of gravitational self-force (GSF) is responsible for the CO's 
gradual inspiral towards the supermassive black- hole. The small parameter n/M allows 
treating the spacetime metric and the CO's orbit perturbatively. At the leading order of this 
approximation the orbit traces a geodesic in the background spacetime of the supermassive 
black- hole Q]. At the next order, the CO's interaction with the 0(/i) first-order metric 
perturbations (FOMP) produces a first-order GSF that accelerates the CO in the background 
spacetime j^, [J. At the next order, the interaction with the 0([i 2 ) second-order metric 
perturbations (SOMP) produces a second-order GSF, and so on. 

Extreme mass-ratio binaries (e.g., Mj /j = 10 5 ) are valuable sources of gravitational waves 
GW) that could be detected by the planned Laser Interferometer Space Antenna (LISA) 
J|. Detection of these sources and determination of their parameters using matched- filtering 
data-analysis techniques requires preparation of gravitational waveform templates for the 
expected GW. Here one of the main challenges is to calculate the GW's accumulating phase, 
e.g. a wave-train from one year of inspiral can contain about 10 5 GW cycles [5lj. Successful 
determination of the binary parameters using matched-filtering techniques often requires 
GW templates in which the phase error is less than one radian over a year 

The following simple scaling argument reveals how the phase scales with /x, and allows us 
to quantify the desired degree of accuracy for the metric perturbation expansion (see also 
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. For simplicity consider a CO that gradually inspirals between two otherwise circular 
geodesic orbits in a strong field region of a Schwarzschild black hole. We are interested in 
estimating the effect of the GSF on the accumulated phase of the emitted GW. Due to the 
GSF the orbital frequency slowly changes from its value at an initial time, and after time t has 



elapsed this frequency has shifted by approximately Cot from its initial value, where Co 



du) 
dt ' 

and t denotes the elapsed time (from initial time) in Schwarzschild coordinates. We denote 
A0 the part of the phase shift of the GW (between two fixed times) which is induced by the 
shift in the orbital frequency. Recalling that the GW frequency is proportional to the orbital 
frequency we find that after an inspiral time Aij ns , the phase shift A0 is approximately 
proportional to Atj ns Co. We shall now find how the quantities in this expression scale with 
ji. Let us consider first Co. We write this quantity as Co = fpE, where E denotes the particle's 
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energy per unit mass. Notice that ^ is independent of \i (it is obtained from an analysis of 
a circular geodesic worldline) while E is determined from the four-acceleration and therefore 
depends on fi. The first-order GSF produces the first-order terms in the expansions of E 
and Co; these terms are denoted here Ei and Co\, respectively Since E\ scales like \i we find 
that loi is 0(/iM~ 3 ). At the leading order the inspiral is driven by the first-order GSF and 
therefore the inspiral time Ati ns is of the order AEE^ 1 = 0(M 2 yU _1 ), where AE is the 
energy difference between initial and final orbits. After Ati ns the term Co\ will give rise to 
a phase shift of order At^ ns Cox = 0(M//i). The second-order GSF produces second-order 
terms in the expansions of E, and Co. These terms are denoted here Ei and C02, respectively. 
Since E2 scales like /i 2 it gives rise to an C02 which is 0(/z 2 M~ 4 ). After Ati ns the term C02 will 
produce a phase shift of order At 2 ns cj 2 = 0[(M/ fi) ]. Therefore, a calculation of A<p to the 
desired accuracy of order (M//z)° (needed for LISA data analysis) requires the calculation 
of the CO's interaction with its own SOMP. 

The challenge of constructing long waveform templates with 0[(M//z)°] accuracy provides 
a practical motivation for the study of SOMP in this article. Furthermore, construction of 
SOMP will extend the applicability of the perturbative scheme to binaries with smaller Mj /1 
mass-ratios. 

Consider the limit where the spatial dimensions of the CO approach zero; we shall refer 
to this limit as the point particle limit (Below we shall make more precise definitions of the 
CO and the limiting process that we use.). In this limit the SOMP away from the CO satisfy 
the following equation 

D[l] = VhVhkhWh. (1) 

Here we used a schematic notation, where the SOMP are denoted /i 2 Z, D denotes a linear 
partial differential operator, the retarded FOMP are denoted V denotes a covariant 
derivative with respect to the background metric, and & denotes "and terms of the form 
In the point particle limit (certain components of) the FOMP in the Lorenz gauge 
diverge as r _1 , where r denotes the spatial distance from the object. Here we encounter a 
problem, since the 0(r _1 ) divergent behavior of h implies that the source term of Eq. 
diverges as r -4 . A naive attempt to construct the standard retarded solution to Eq. ((TJ) by 
imposing Lorenz gauge conditions on I, and then integrating the resultant wave equation 
using the retarded Green's function, results in an integral that diverges at every point in 
spacetime. A similar problem in a scalar toy-model has recently been studied |9[. In this 
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article we develop a regularization method for the construction of well-defined and physically 
meaningful solutions to Eq. (JTJ). 

We should mention here another regularization problem coming from infinity. For an 
infinitely long worldline, the leading asymptotic behavior of the source term in Eq. (JTJ is 
0(r~ 2 ). This behavior at infinity renders the retarded solution I divergent 1 . This issue, 
however, lies outside the scope of this article. Here we assume that a regularization at 
infinity had been carried out. 

For simplicity we choose the CO to be a Schwarzschild black hole, which moves in a 
vacuum background geometry characterized by length scales that are much larger than /i. 
Designating the length scales associated with the Riemann curvature tensor of the back- 
ground geometry with {Hi}, and denoting TZ = min{7Zi}, we express the above restriction 
as n <C TZ. 

Our analysis is based on the method of matched asymptotic expansions (see e.g. jlO, 11|), 



in which different approximate solutions to Einstein's field equations are obtained in differ- 
ent, overlapping, regions of spacetime. These solutions are then matched in their common 
region of validity. Consider the following decomposition of spacetime into an internal zone 
that lies within a worldtube which surrounds the black hole and extends out to r = rj(7Z) 
such that Ti <C 7Z, and an external zone that lies outside another worldtube at r = te{h), 
such that jj, <C r^. The interior of this inner worldtube is denoted S. By virtue of fi <C TZ 
we choose rg to be much smaller than 77 such that there is an overlap between the above 
mentioned regions in te < r < rj. We refer to this overlap region as the buffer zone (in the 
buffer zone r can be of order y/JjjZ). 

In the external zone we decompose the full spacetime metric gj^f 1 into a background 
metric (e.g. a spacetime of a supermassive black hole), and a sequence of perturbations 
that are produced by the Schwarzschild black hole (with mass /i), reading 

9Sf(x) = g^{x) + mV(z) + + 0(/i 3 ) . (2) 



1 Asymptoticly h has a form of a gravitational wave. Therefore the dominant terms in an asymptotic 
expansion of V/i decay like r _1 . This implies that the source term in Eq. is asymptotically 0(r~ 2 ). 
Notice that this source term has a static 0(r~ 2 ) component which does not vanish after averaging over 
time. An attempt to construct the retarded solution to Eq. Q with this static term present yields a 
divergent expression (even if we resolve the difficulty with the singularity in the vicinity of the worldline) . 
Ori has recently suggested a resolution to this problem . 
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Here the dependence on /x is only through the explicit powers fi l . Throughout we use 
the background metric g^ v to raise and lower tensor indices, and to evaluate covariant 
derivatives. Substituting expansion (J2J) into Einstein's field equations in vacuum yields 
linear perturbation equations for h^ u and l^ u . These equations are valid for x ^ S. Here we 
shall be interested in the limit \i — > 0. Note that this limit has no effect on the forms of the 
perturbations equations. Yet their region of validity is affected since we let te{^) — > as 
/j — ► 0. In this limit the perturbation equations read 

D lu ,[h]=Q, x$z(t), (3) 
D^] = S, u [h] , x?z(t). (4) 

Here z(r) is a timelike worldline, where r denotes proper time with respect to the background 
metric. An overbar denotes the trace-reversal operator h^ v = h^ u — (l/2)g fMU h a a ; for brevity 
we have omitted tensorial indices inside the squared brackets. D^ u and are obtained 
from an expansion of the full Ricci tensor. By substituting g^u + Sg^ into the Ricci tensor we 
obtain Rff 1 = R { $ + R$±6g] + R$[6g] + 0(6 g 3 ). (For explicit expressions for the terms m 



this expansion see e.g. [12j . We adopt the sign convention of this reference for the Riemann 
tensor). To simplify the notation we denote R$ [h] by D^ iv \h], and denote —Rj$ [h] by S^h], 
where h^ u is expressed in terms of h^. 

By imposing the Lorenz gauge conditions h^ u . u = 0, and matching with the internal-zone 
solution, D'Eath showed |l| that h^ v is identical to the retarded FOMP produced by a unit- 
mass point particle tracing the worldline z(r). At this leading order of approximation z(r) 
is a geodesic of the background spacetime 1|, denoted zg(t). Substituting z(r) = zq(t) in 
Eq. (J2J) provides us with sufficient accuracy for the calculation of the FOMP. However, to 
calculate the SOMP we must account for the 0(/x) acceleration of z{r) due to the first-order 
GSF corrections. To account for these corrections we use the gauge dependence of the first- 
order GSF Ijjj], and choose a convenient gauge in which the first-order GSF vanishes. In 
this gauge, the geodesic worldline zcij) is sufficiently accurate for the construction of the 
SOMP. 

To spell out the desired gauge conditions we examine the expression for the 0(/i) accel- 
eration which is produced by the first-order GSF in the Lorenz gauge ^| 

^ = -f(<r + uVVti'pVp/i* - V v h%) . (5) 
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Here all quantities are evaluated on the worldline. We use the Detweiler- Whiting decom- 
position hfj, u = h^ u + h^ u , where h^ u and h^ u are certain regular and singular potentials, 
respectively (for their definitions and properties see Consider now a regular gauge 

transformation generated by a vector field In the new gauge the FOMP, h^ u , are given 
by = h % + hfy , where h% F) = hf v + + £„ ;jtl (The gauge invariance of h^ v follows 
from the analysis in Ref. Q].). 

Many gauge choices can provide us with the desired requirement d 1 = 0. For simplicity 
we impose the following gauge conditions on the worldline: 

[*2"U, = °. WW- (6) 

We shall refer to this gauge as the Fermi gauge. We may now replace the restrictions x ^ z(r) 
with x z G {r) in Eqs. (|3I4|) . Notice that the gauge conditions (JHJ) are specified only on 
Zg{t). To complete the gauge construction we adopt an arbitrary continuation of £ M to the 
entire spacetime (e.g., continuation along future null-cones based on zg{t) in a manner that 
preserves causality 2 ). 

Before solving Eq. (JIJ) we study the singular properties of near zq (r). Consider 
expanding in the vicinity of zg{t), on a family of hypersurfaces r = const that are 
generated by geodesies that are normal to the worldline. We use Fermi normal coordinates 
based on Zg(t), thus obtaining simple expressions. We find that 

S, v (x) = A, u r- 4 + 0(r- 2 ). (7) 

Here = denotes equality in Fermi normal coordinates. We have A^ u = Au^u u + Irj^ — 
14Q^Q U , r}fj, v denotes the Minkowski metric, r = \f5 a bX a x b , where x a denote the spatial 
Fermi coordinates (a, b take the values {1,2,3}); — 5q is a vector field, fi a = x a /r, 
Cl° — 0, and VL V — g vp VL^. The absence of 0(r~ 3 ) terms in Eq. (JJJ) follows from the gauge 
conditions © together with the fact that zg(t) is a geodesic worldline. 

2 To construct this continuation consider the future null-cones E T emanating from the worldline Zg(t). 
Here we focus on a local neighborhood of the worldline in which these null-cones do not intersect each 
other. For an arbitrary point zq{t~) one may choose an arbitrary continuation of £ p on E T -, such that 
£^ decays to zero away from the worldline. In this way the constructed gauge preserve causality in the 
following sense: The perturbations on the null-cones S T for r < t~ will remain unchanged if one modifies 
the worldline for r > t~ . Such a modification of the worldline is possible by introducing additional GW 
that interact with the worldline for r > t~ . 
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To tackle the r 4 singularity of we decompose into two tensor potentials, reading 

Itxv = Vv + > ( 8 ) 
where satisfies D^ip] — A^ v r~ 4 + 0(r~ 2 ). Notice that this equation has the same r~ 4 
singular source term as Eq. P|). Here, however, we do not impose restrictions on the terms 
of order 0(r~ 2 ). Suppose that we can construct a solution ip^, then by subtracting D^ty] 
from both sides of Eq. (jlj) we obtain 

D^[«l = ^ ,x£z G (r). (9) 

Here 5S^ U = S^ u — D^[ip}. By construction diverges as r -2 , while S 1 ^ in Eq. (JJJ) 
diverges as r~ 4 . In this sense Eq. © is simpler than Eq. (|1J). 

To construct f/y^ we use a linear combination of terms that are quadratic in h F . Since 
h F u diverges as r _1 we find that by applying the operator D^ v to terms which are quadratic 
in h F u we obtain terms which diverge as r -4 . First we construct four independent quadratic 
tensor fields: ^ = h F %h^, < = h F %h F v , ^ = (h F ^h F Jg^, tf u = (h Fp p ) 2 g^. 
Combining these terms we may construct the desired ip^, which reads 

= l~ A [2(ca< + c B <f%,) - 7(c c tf v + c D <f%,)] • (10) 

Here the constants ca, cb, cq, cd must satisfy ca + cb = 1, cc + = 1, but are otherwise 
arbitrary. Roughly speaking, ip^ captures the asymptotic behavior of the Schwarzschild 
solution at second order. Due to this property ip^ is an exact solution to Eq. (J3J for the 
case of a flat background spacetime. 

Having constructed ip^, we now face the problem of solving Eq. Q. For this purpose we 
invoke a purely second-order gauge transformation of the form x^ — > x p — /i 2 £|^ that allows 
us to impose the Lorenz gauge conditions 5l fJ,u . u = 0. First we seek a particular retarded 
solution to Eq. (j^J). Here it is useful to remove the restriction x zq(t), and continue 
SSpv to the worldline zg(t). Clearly a solution to the continued equation also satisfies the 
original equation (i.e. with the worldline excluded). We choose the simplest continuation 
by demanding that no additional singularities (e.g. delta functions) are introduced on the 
worldline i.e., the only singularities on the worldline must be those originating from a local 
expansion of S^ u . Eq. (jHJ) now reads 

USl^ + 2Ry v 8l vp = -25S,„ . (11) 
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Here □ = g aP V a Vp. We define 81^ to be the retarded solution to Eq. (fTTj). Since 8S^ U 
diverges only as r~ 2 its retarded solution has a finite contribution originating from the 
vicinity of z G {j). 

By definition the retarded solution 81 ^ u satisfies Eq. (111)) . However, we still have to 
verify that it also satisfies Eq. Q. This equation will be satisfied if 81 ^ satisfies the Lorenz 
gauge conditions. To investigate this point one can apply the divergence operator to Eq. 
fITTj) . thereby constructing a differential equation for VSl^. The fact that VSl^ = then 
follows from the properties of the source term of this equation, namely — 2V V 5S' M!/ . Using a 
perturbation expansion of the Bianchi identities one finds that VSS^ = for x ^ zg(t). 
This property together with an analysis of VSS^ as one approaches the worldline reveals 
that the retarded solution to Eq. (fTTj) satisfies the Lorenz gauge conditions, and therefore 
it also satisfies Eq. Q as required. (The full analysis will be given elsewhere 

So far we have constructed a particular solution (JHJ) to Eq. We now construct the 
general solution 1^ to this equation. Later we shall impose a set of additional requirements 
on ljj vJ and thereby obtain a particular solution which is physically meaningful. Since Eq. 

is valid for x ^ zg{t), we find that we can construct a new solution by adding to 
a potential that satisfies a semi-homogeneous equation ,i.e., a homogeneous equation for 
x ^ zg(t), reading 

D flu [l SH ] = ,x#zg(t). (12) 
The general solution to Eq. (JH) is given by = I^f + T„ v , where l s J^ is the general solution 
to Eq. (EJ). 

The set of additional requirements to be imposed on I® can be expressed as requirements 
on l s J^ . Denoting 7^ the semi-homogeneous potential that satisfies these requirements, we 
express the desired physical solution W v as 

^ = l^ + 7^ = $iw + M l u, + 7i»>- (13) 

We consider the following additional requirements: (i) gauge conditions, (ii) causality re- 
quirements, (iii) global boundary conditions, and (iv) boundary conditions as x — > z g {t). 

(i) We impose Lorenz gauge conditions on 7^. Eq. (|T2*|) now takes the form 

n% v + 2R\\% p = , x £ z g {t) . (14) 

(ii) We define 7^ to be a retarded solution to Eq. (JHJ). Recall that 8L U is the retarded 
solution to Eq. (fTTj). and ?/> M „ is completely determined by hF v . This construction implies 
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that if we prescribe initial data for the metric perturbations on a spacelike hypersurface So, 
then the physical solution l^ u (x) at the future of So is unaffected by an arbitrary modification 
of the initial data outside J~(x) fl S , and in this sense the second-order solution tf u (x) is 
manifestly causal, (iii) We demand that the only source of 7^ is the black-hole itself, thus 
excluding any incoming waves, (iv) Let us consider once more a small but finite mass //. 
Here, Eq. (J14|) is valid at x $ S. Following D'Eath's analysis of FOMP we express a 
solution to Eq. (jl4j) using a Kirchhoff representation. Denoting S# the boundary of S, we 
express 7^ in terms of the following surface integral 

r [ feWM^ ve V'^V) (is) 

-^P\x f )V e 'G^ a , p ,[x\x f ])dX e ,. 



Here G T f in i \x\x'] denotes the retarded Green's function (for its definition and properties 



see e.g. ll|), c?S e / denotes an outward directed three-surface element on Consider 
substituting an expansion (in powers of r) of a given potential 7 Q/3 into Eq. (|15|) . and then 
taking the limit fi — > 0. Recall that in this limit we have te{h) — > 0, and notice that G?S e 
scales as r\. Therefore, only the diverging terms (as r — > 0) in the expansion can contribute 
to ^i^ v {x) in the limit. These divergent boundary conditions for Eq. (fT4"]) are obtained from 
Eq. (fT3j) together with an analysis of the divergent behavior of Tp^, 51^, and in near 
r = 0. 

First we expand ip^ near r = and obtain ip^ v — —-^[2u^u v + 77]^] + O(r ). Next we 
consider 51^. Solving Eq. (|TT|) approximately in the vicinity of zg(t) reveals that 51^ is 
bounded as r — > |l6|. Finally, we expand using the internal-zone solution, which is 
discussed next. 

We now calculate the desired diverging terms in the expansion of 1^. For this purpose 
we use the method of matched asymptotic expansions. For simplicity we assume that all 
the length scales characterizing the background spacetime {Hi} are of the same order of 
magnitude TZ. Expanding the full spacetime metric in the internal-zone using the smallness 
of r/TZ and fJ.fR. yields 



& + rc-^g? + ir 2 g<$ + o(n- 3 ) . (16) 



Here gf^f is the metric of the Schwarzschild black-hole. Recall that in the buffer-zone both 
Eq. (Jrfij) and Eq. (J2J) are valid, and should therefore match with each other. To match 



the internal-zone metric with the external-zone metric we further expand Eqs. (|2I16J) in the 
buffer-one. 

First we consider the internal-zone expansion. The first term in Eq. is the 

Schwarzschild metric g^ h . In isotropic cartesian coordinates this metric takes the form 



ds 2 = -{2p - p) 2 (2p + fi)- 2 dt 2 + [1 + p(2p)- 1 } 4 (dx 2 + dy 2 + dz 2 ), where p 2 = x 2 + y 2 + z 2 . 



In the buffer-zone p <C p and therefore the Schwarzschild metric can be approximated with 
an asymptotic expansion in powers of pp~ x giving 



where H 1 = 2(77^ + 25j5°), H 2 = 1/2(3^ - 6° 6$. The next term in Eq. JEJ is TZ^g^J. 



In a suitable gauge this term vanishes identically [10]. 

Next we consider the external- zone expansion Recall that in the buffer- zone r<K. 
We therefore expand the terms in Eq. (J2J) in powers of rTZ~ l . The first term in Eq. (J2J) is 
the background metric g^. Expressing this metric with Fermi normal coordinates based on 
the worldline, and expanding it in powers of r1Z~ l gives g^ — i]^ + 0(lZ~ 2 r 2 ). Expanding 
the next term in (J2J gives ph^ — 2pr~ 1 (r]^ u + 2u IM u u ) + 0(p1Z~ 2 r 1 ), and expanding the third 
term p 2 l^ u gives a sum of a schematic form p 2 5^°^ 7Z~ %r% ~ 2 ■ Recall that here we are only 
concerned with divergent boundary conditions that are encoded in the % — 0, 1 terms. 

First let us consider matching the first three terms in the expansion of the Schwarzschild 
metric 1)170. all these terms scale like TZ°. From the above expansions we find that both g^ h 
and g^y have the same leading term - rj^. The next term in Eq. ()17)) is /ip _1 if this term 
coincides with the leading order expansion of ph^ u , where we identify p with r. The next 
term in Eq. (|T?j) provide us with the desired i = term in the expansion of p 2 lj^ u - We find 
that this term is identical to the leading term in the expansion of ip^. 

Next let us consider matching the terms that scale like 1Z~ X . Since 7Z' 1 g$ = all the 
terms that scale like 1Z~ 1 vanish. This conforms with the fact that the expansions of g^ v 
and ph^ u do not contain terms that scale like Tlr 1 . This also implies that in the expansion 
of p 2 l^ v the term with the schematic form p 2 TZ~ 1 r~ 1 is zero. Combining the i = and % = 1 
terms in the expansion of p 2 l^ u we conclude that — —^2[2u^u u + 7r]^ v ] + 0(r°). 

Combining the expansions of and near r = 0, together with the fact that 51^ 
is bounded in the vicinity of the worldline, and using Eq. (|T3j) we find that 7^ does not 
diverge as x — > zq{t). As was previously discussed only divergent boundary conditions (as 



9, 



: Sch 



% v + pp 1 H] lu + p 2 p 2 Hl u + 0(p 3 p 3 ), 



(17) 
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r — > 0) can produce non- vanishing contributions to 7 MV (x). Since we found that 7^ does 
not diverge in the vicinity of the worldline we conclude that 7 M1/ = 0. By virtue of Eq. (JT3Jl 
we finally conclude that the physical SOMP are given by 

jJL = VV + ^ • (18) 

Here ip^ u is given by Eq. (JTUJ), and 51^ is the retarded solution to Eq. (JTTJ). By construction 
satisfies the previously mentioned requirements. In particular it is manifestly causal, 
and it matches with the internal-zone solution. Moreover, Eq. (|18j) provides a simple 
covariant prescription for the construction of the SOMP, without any reference to a particular 
(background) coordinate system. This allows a considerable amount of flexibility in the 
construction of the SOMP. 
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